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ABSTRACT 

We use the canonical formalism developed together with David Robinson to study the Einstein equations 
on a null surface. Coordinate and gauge conditions are introduced to fix the triad and the coordinates 
on the null surface. Together with the previously found constraints, these form a sufficient number of 
second class constraints so that the phase space is reduced to one pair of canonically conjugate variables: 
A 3 2 and S 3 2 . The formalism is related to both the Bondi-Sachs and the Newman-Penrose methods of 
studying the gravitational field at null infinity. Asymptotic solutions in the vicinity of null infinity which 
exclude logarithmic behavior require the connection to fall off like 1/r 3 after the Minkowski limit. This, 
of course, gives the previous results of Bondi-Sachs and Newman-Penrose. Introducing terms which fall off 

■ more slowly leads to logarithmic behavior which leaves null infinity intact, allows for meaningful gravitational 
radiation, but the peeling theorem does not extend to in the terminology of Newman-Penrose. The 

■ conclusions are in agreement with those of Chrusciel, MacCallum, and Singleton. This work was begun as a 
preliminary study of a reduced phase space for quantization of general relativity. 

^ ' 1 T / 

. 1. Introduction. 

The canonical approach to quantum gravity received a strong impetus from the introduction of the new 
variables by Abhay Ashtekar [1-4] . The use of a self-dual connection one- form and a vector density triad as 
canonical variables leads to a Hamiltonian which is a polynomial of degree four in the new variables. This 
structure suggests a simplification of the canonical formalism which may lead to significant easing of the 

l/"") | problems of quantization. By now there is a considerable literature detailing the efforts of many people to 

OS ■ understand how to make effective use of these variables [3,4]. 

""q I As part of this overall effort, together with David Robinson, we have renewed the study of canonical 

general relativity on a null surface [5]. This program had previously been undertaken using the metric or 
5^ | a tetrad as the configuration space variables [6,7]. Although these efforts did not recover all the Einstein 

equations in a natural way, their principal drawback is that the resulting system of second class constraints is 
complicated by the non-polynomial structure so that there does not appear to be any hope that a successful 
Dirac quantization [8,9] can be carried out. Our hope was that use of the self-dual formalism on a null 
surface would retain the polynomial structure of the Hamiltonian and the constraints. If so, the second class 

■ constraints might not be as formidable as in the previous treatments. 

While some of the second class constraints are indeed simpler, for the most part they are sufficiently 
complicated that they cannot easily be eliminated either directly or by means of the introduction of Dirac 
brackets [8-10]. Therefore, as a first step toward considering a reduced phase space quantization, we have 
repeated the analysis of the gravitational field in the vicinity of null infinity. The use of the Ashtekar variables 
for this analysis falls between the metric formulation of Bondi [11] and Sachs [12] and the spin coefficient 
method of Newman-Penrose [13,14]. Although the full complement of Dirac brackets cannot be obtained, a 
machine caculation did show that the bracket of one component of the connection (related to the expansion 
and shear) with its conjugate triad density is the same as the corresponding Poisson bracket. This identified 
for us which variables should serve as the true dynamical variables in terms of which to express the remaining 
components of the gravitational field. 

By focussing attention on the connection, we were immediately struck by how deep in a 1/r expansion 
one must go in order to avoid the appearance of logarithmic behavior. Therefore, the question naturally arises 
whether such logarithmic behavior is consistent with asymptotic flatness and finite energy. While we were 
still analyzing this question, Chrusciel, MacCallum, and Singleton [15] showed that this is indeed the case. 
However, because in this work we use a different formalism in which the need to examine the logarithmic 
terms is glaring, it is worthwhile to present the results. Also, there is a difference in our solutions. We find 
that even with the appearance of logarithmic behavior, the coefficient of 1/r in the connection must vanish 
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whereas they do not. A detailed discussion of our solution leading to this difference is left for the Appendix. 

To see how far the connection must drop in 1/r in order to avoid logarithmic behavior, we shall first 
carry out the analysis requiring an expansion in 1/r without logarithms and then go back to indicate what 
changes would be introduced by the logarithmic terms. The simplification in the calculations by neglect of 
logarithms is enormous. In agreement with the earlier work of Novak and Goldberg [16, 17], the present 
results and those of Chrusciel, MacCallum, and Singleton show that null infinity can be defined and that 
energy-momentum and radiation of energy-momentum remain finite. In our case, logarithmic terms come in 
below those responsible for radiation. In the more general considerations of Chrusciel et al, the coefficients 
of the leading logarithmic terms are independent of the time. This point does not come up in our work. 
However, in both studies, the logarithmic terms come in before those needed to define angular momentum. 

In the following section, we shall give a streamlined review of the construction of the Hamiltonian on a 
null surface [5]. In section 3, we introduce our coordinate and tetrad conditions and present our analysis of 
the equations, listing the order in which they are to be solved. The solution in the absence of logarithmic 
terms is given in section 4 and in section 5 we discuss the logarithmic terms. We close with a discussion of 
our results. 
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Section 2. The Hamiltonian. 

To obtain the Hamiltonian, we started from a complex Lagrangian constructed from the self-dual part of 
the Riemann tensor, following similar work by Jacobson and Smolin [18] and by Samuel [19] on a space-like 
surface. We assume the space-time to be real, but consider complex solutions of the Einstein equations. 
After the calculation has been completed, we impose reality conditions on the variables and recover a real 
metric and curvature. Because a null surface is degenerate, we lose one of the field equations if we allow 
the initial surface to be null from the beginning. Therefore, in order to be sure that we recover all of the 
Einstein equations, we included an auxiliary variable a. The surfaces t =constant are space-like, time-like, 
or null, when a < 0, > 0, or = 0, respectively. To guarantee that the surfaces would be null, we adjoined 
a 2 = to the Lagrangian with a Lagrange multiplier. Moreover, because working on a null surface imposes 
constraints not present on a space-like surface, we did not a priori eliminate any variables as non-dynamical. 
Therefore, we started with a phase space of 40 variables. The 12 first and 14 second class constraints leave two 
dynamical degrees of freedom per hypersurface point as is appropriate on a null surface [20] . An examination 
of the resulting structure shows that we may limit the phase space variables to the nine components of the 
connection and of the densitized triad vectors on the hypersurface. We shall do so below. 

We introduce the null basis of one forms and the dual tetrad basis 

6° = (N + av 1 l N l )dt + av 1 l dx i , 

9 i = t/^dt + dx*) 

e = N-^dt-Wdi), 

ei = {vi i + aiN- 1 N i )d i -aiN- 1 dt, 

where cti = a5 1 - 1 and v l ivf — (5'j. All indices have the range 1-3 and repeated indices sum. Bold face 
indices refer to the one forms and tetrads. The signature is -2, with 

ds 2 = e° ® e 1 + e 1 ® 6>° - e 2 ® e 3 - e 3 ® e 2 , 

= Va0O a ® 

(Greek letters range and sum from zero to three.) It follows that the surfaces t ^constant are null surfaces 
when a = 0. 

We choose the orientation so that the volume form is 

V = -i6° A9 1 A0 2 A9 3 . (2.2) 

Thus, —iNv is positive, where v is the determinent of i/j and the components of the Levi-Civita tensor with 
respect to this basis, e Q( 3 7< 5 and e al3jS , satisfy 

(2.3) 

The connection coefficients are defined by 

d0 a = 

and the Riemann tensor by 

dw a p + w" 7 A u~< [} - ^R a /3-ysO J A 8 s =: R a p 
We take as a basis of self-dual two-forms 

S 1 = ^(e 1 A6°+6 3 A6 2 ), 

S 2 =9 1 A6 2 , ( 2 - 6 ) 
S 3 = e 3 A6°. 
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(2.1) 



(2.4) 



(2.5) 



These two-forms satisfy the equations (A, B = 1,2, 3) 



*S A = iS A 
S A f\S B = ig AB V. 



The S*0(3) invariant metric 



and its inverse, 
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are used to raise and lower the uppercase Latin self-dual, triad indices. 
The self-dual components of the connection are 

and are represented by T A 7 : 



(2.7) 



(2.8a) 



g AB = -1 (2.86) 



T 1 ^ i(^ 01 +^ 23 ), T 2 :=uj 21 , T 3 :=LU 03 . (2.9) 
From these we obtain the self-dual components of the Riemann tensor as 

^R A = dT a + ri A B C T B AT C . (2.10) 
In a 3 + 1 decomposition, we have 

T A = A A l dx'+B A dt, (2.11) 

and 



(2.12) 



R A l3 = 2A A {iJ] + 2 V A BC A B j A c i , 
R A oi = DiB A — A A ifi. 
The derivative operator D acts on the index A as 

DJ A := 0J A + 2 V A BC A B J C . 

The complex action introduced by Jacobson and Smolin [18] and by Samuel [19] is 

S = J R A AS B g AB . (2.13) 

With the above decomposition, the Lagrangian density formed using the self-dual curvature tensor takes the 
form 

C = (A A ifi E A i + B A D^ A * + R A t3 N^ A i - 

- Nv\R\ 3 ^ + fl 2 «£i J ')) + ^(ly - <«>') + P(uf- (2 ' W) 

In the above we have introduced v l := f5 2 ' and 

Si* := -vj*, S 2 J := -ai 2 \ £ 3 * := -tfe*. (2.15) 
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The over (under) tilde indicates a density of weight plus (minus) one. Thus, v l and T,a 1 are densities of 
weight one while the other variables are not dcnsitized except as indicated by the tilde. 

It is evident from this form of the Lagrangian that only variation of A A i and Y<a 1 lead to dynamical 
equations. The remaining variables lead to equations on the hypersurface t = constant, hence to constraint 
equations. In the absence of coordinate or tetrad conditions, TV* and B A are arbitrary while a and ^ are 
needed to make t = constant a null surface and to obtain the otherwise missing field equation. The strange 
thing is that Nv l appear in that combination and cannot be solved for separately. As long as the solution 

is complex, N is arbitrary. But, when the reality conditions are imposed, v l = £ 3 l and N will no longer be 

undetermined. This structure in the Lagrangian foreshadows the fact that the scalar constraint is second 
class. 

Therefore, we take the phase space to be defined by (A A i7 S^ 1 ). It is easy to see that they form 
canonically conjugate pairs with the Poisson brackets 

{A A i {x), EbV)} = S A B 5 3 (x-x'). (2.16) 

Then the Hamiltonian takes the form 

H = J d 3 x{NH + N*H t - B A Q A - - p(a) 2 } (2.17) 

where 

Ho :=v l {R\^ +R 2 ij ^) = 0, 



Hi 
Qa 



-R A i^ A j = 0, 



(2.18a) 



n 

D^ A l = 
E 2 * + m/ = 0. 

are constraints which arise from varying the action with respect to N, N l , B A , and 
Hi. respectively. Other constraints come from varying p, a, and v l : 

a = 0, v% = 0, (2.186) 

and 

(j> t := + = 0. (2.18c) 

Propagation of the constraint C l = leads to 

X*:= 2S b 2 D j ^Nv^a j] Q A b) -2A 3 J N^ 1 j ^ - B 3 ^ = 0, {2A8d) 
Q A b := SiS^ + S A S%. 

These are all the constraints, but propagation of the constraint Q 3 ~ leads to a condition on /Ltj which is 
important in obtaining the otherwise missing equation: 

- R\ ^3 - R 1 ^^. (2.19) 

Note that v 1 ^ = H and SiVi = ^H-i- 

Thus, there are 14 constraints among the dynamical variables and three conditions on the Lagrange 
multipliers a and /Ltj. Five of the constraints, Hi,Gi, and Q2 are first class while the remaining constraints, 
including Ho are second class. Three of the second class constraints are conditions on v l , so there are 16 
conditions on the 18 phase space variables per hypersurface point. 

The Hamiltonian equations of motion for the dynamical variables are 

A x h0 = 5 1 A DiB A + N'R 1 i j -Nv i I^i j , (2.20a) 
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A\ = 5 2 A D t B A +NiR 2 ij- m, (2.206) 
A 3 ifi = S^aD^ + N^R^j-Nv'R 1 ^, (2.20c) 

Ei^o = 2S B 1 D 3 (Nv^ i Y; A rtQ A B ^ -25 B 1 D j (N^ B A) -2B 3 Z 3 i , (2.20d) 

E 3 \o - 2<5 B 3 ^ (iV^E^]Q A B ) - 2<5 s 3 ^(iV[ J E B ^) 

V ~ ' (2.20e) 

+ 2 J B 1 E 3 J + J B 2 Ei\ 

The quantities we have introduced are complex. To recover the Einstein theory we must impose reality 
conditions at some point. These conditions are that 

Ei* = u* = E 3 \ 

and that 

r 1 -2co 23 , F 2 =u; 31 , andf 3 ^ 02 
form the anti-self-dual components of the connection. 
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3. Analysis of the Equations. 

We assume that space-time is asymptotically Minkowskian and that outside of a timelike cylinder, the 
coordinate t defines a congruence of hypersurfaces. When a = 0, these are null surfaces which have the 
topology of null cones extending to null infinity. These null cones in turn are foliated by closed two-surfaces 
so that each null generator is labeled by the usual angular coordinates (6, (f>) of the unit sphere. Following 
Bondi and Sachs [11,12], we choose the coordinate r along the generators to be the luminosity distance; 
that is, the area of each two surface r = constant is Airr 2 . The coordinates x l of the previous section are 
then (r,0,(j>). For the boundary conditions, we assume that T,a 1 and v l , thus the metric as well, take their 
Minkowski space behavior in limit of null infinity. In as far as it is possible to do so without losing the 
gravitational radiation, the same is true for the connection. It turns out that A 2 i contains radiative terms 
in this limit. These conditions, including the radiative terms, are unaffected by super-translations. 

The five first class constraints allow us to make a convenient choice for the triad densities and the 
coordinates. One coordinate condition has been used to fix r as the luminosity distance: 

-ir 2 sin9 r] ljk v j T, 3 k = v 2 . (3.1a) 

With the other two available coordinate conditions, we can set 

Si 1 = -ir 2 sinO 5\. (3.16) 

Then the four real functions in the null rotations generated by Q\ and Qi allow us to fix S 3 * tangent to the 
surfaces r =constant, that is, S3 1 = 0, and then to set A x \ = 0. The latter is equivalent to setting e = in 
the Newman-Penrose formalism [14]. 

As pointed at the end of the previous section, the Dirac brackets of the conjugate pair (A 3 2 , S3 2 ) are 
equal to their Poisson brackets. Therefore, these can be identified as the dynamically independent degrees of 
freedom for the gravitational field. However, they are not observables because they arc not diffeomorphism 
invariant. Nonetheless, they represent the initial data which can be specified on an initial null surface. The 
remaining variables and parameters are then determined by any remaining gauge freedom, the constraints, 
and relations from the propagation equations. 

With the coordinate and gauge conditions given above, the constraints and propagation equations have 
a natural order for their solution. Below we will give the equations in the order in which they can be solved. 
Which variable is to be solved for is generally clear. (The indices a, 6, .. range over 2 and 3.) 

Gx = 0^ 1 \ 1 +2A 3 l Z 3 t = (3.2a) 

£2 = =^A 3 jEi* = (3.26) 

U x =0^A\ 1 £ 3 J = (3.2c) 

Q 3 = ^£ 3 \i - 2A 1 l S 3 i - A 2 jSi* = (3.2dl) 

7~ia = =>\A^~ a .j — A^j^a — A 2 a A 3 j + ^4 3 a ^4 2 j]5]l- , 

+ [A 3 a , 3 -A a jia + A a a A 1 j -A a j A 1 a ]^ - 

[A^j - A 1 j ta - A 2 a A 3 j ]Y, 1 j + A a a \S a \i - 2A\^] + 

[A a a>j -A a j!a + A a a A 1 j -A a j A\]^ = (3.2d2) 

a = ^A 2 aj 5V - Z 3 c A 3 c A 2 b Z b a + G a = (3.2e) 

V b A 3 

G a := (A 1 ^ - A\, a )E 3 b - {A 2 j, a + 2A 2 ,A\)^ 
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x a = =*(V IV) j + Z 3 c A 3 c Z a b V b = 0, (3.2/) 
Nv a = V a 

X 1 ={)^B 3 Y, 1 1 + {V b Ys 1 1 ) :b + 2A 1 b V b Y, 1 1 - A 3 ^^ 1 = (3.2. 9 1) 

± ± a = ^(A^IV'Xj - (T/ a E 3 b - V b Ys 3 a ),b + 2A 2 j Z 1 j 

+ 2A 3 b Z 3 b N a + 2B 3 E 3 a = (3.2.g2) 

Ei 1 = 0^2N 1 A 3 b X 3 b -A 2 j V j 'Z 1 1 -(N a i: 1 1 ), a = (3.2h) 

A 1 ! = ^S 1 ,! + A 2 ii? 3 - A ra (A 1 a: i + A 3 a A 2 i) 

-V a (A\^ AKl -2A\A\) = (3.2i) 

E3 1 = ^B 2 ^ 1 - (A^IVXa + 2iV 1 J 4 1 a E 3 a 

-iVM 2 ^ 1 = (3.2j) 

This completes the set of equations which can be solved cither by integration along the null generators of 
the surface t =constant or algebraically on that surface. The equation forA 3 i = yields another equation 
for B 3 which is then trivially satisfied. 

The remaining equations all contain time derivatives. 

A\ =B\ a + A 2 a B 3 - A 3 a B 2 + N j [A 1 a j - A 1 ^ a + A 2 ^ 3 a - A 3 ^ 2 a ] 

- V b [A 2 a . b - A 2 b . a + A\A 2 a - A 2 b A\] (3.3a) 

A 2 , = B 2 , l -2A 1 l B 2 + 2A 2 l B 1 

+ W [A 2 ij - A 2 js - 2A 1 j A 2 i + 2A 2 jA 1 i ] - m (3.36) 

A 3 a = S 3 , a - 2A 3 a B 1 + 2A\B 3 + N j [A 3 aJ - A 3 jiU + 2A 1 j A 3 a - 2A 3 1 A 1 a ] 

- V b [A 1 b, a - A\, b + A 2 a A 3 b - A 3 a A 2 b ] (3.3c) 

t 3 a = - 2<9 7 -(7V[°5V']) + 4A 1 j N^ a E 3 ^ + 2A 2 j N^Y, 1 A 

-2B 1 S 3 a + B 2 S 1 a (3.3d) 

The fact that the Poisson brackets of the constraints with the Hamiltonian vanish modulo the constraints 
themselves tells us that the propagation of the variables is determined by the propagation of (A 3 2 , S3 2 ) 
alone. That is not quite true because, as in the case of Bondi-Sachs [11,12] and Newman-Penrose [13], 
integration along the null generators introduces arbitrary functions of (t, 0, <p). We find that propagation of 
A 1 a and A 2 a v a lead to the conservation equations for angular momentum and mass. The time derivative 
of A 2 Ei° defines /z^Ei 8 which when equated to (2.19) yields an identically satisfied field equation. On the 
other hand, the time derivative of A 2 a Y> 3 a defines [iiT, 3 l as the null component of the conformal tensor. This 
is in complete agreement with the previous work. 

The integration along the null generators can be carried out from the time-like cylinder to null infinity 
without an expansion in 1/r in a manner similar to that of Tamburino and Winicour [21], but this formal 
result does not exhibit the presence or lack of logarithmic behavior. Therefore, in the next section we shall 
set up the calculation of the asymptotic behavior and then list the results. 
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4. The Asymptotic Solution. 



In this section we shall first solve the complex equations for the triad and the self-dual connection. Then 
we shall apply the reality conditions and show the explicit relationship of our results to those of Sachs [12]. 
We shall look for solutions for which the triad differs from its Minkowski space value by factors with an 
expansion in 1/r. Although there are logarithmic terms consistent with the assumption of asymptotically 
Minkowskian behavior, in this section we shall choose the powers of 1/r to avoid their occurrence. The 
appearance of terms in (In r ) m /r n will be discussed in the following section. The Minkowski space solutions 
are given in Appendix 1. 
The Solution. 

As noted previously, the initial data on the surface t = will be given by (A 3 2 , S3 2 ). By our coordinate 
conditions we have 

Si* = -ir 2 sin 6 8\. (4.1) 

Then, g x and 2 , Eqs.(3.2a,b) tell us that 
To proceed, we write 

A\ 

and find from (4.2a) and Hi that 

A 3 

a 

^3 

The forms (a a , b a ) and the vectors (s a ,t a ) are the eigen- forms and eigenvectors of the flat space part of the 
matrix Z a \, as defined in Appendix 2. They satisfy the following algebraic relations: 

a a s a = b a t a = 1, a a t a = b a s a = 0. 

We assign a spin- weight of -1 to b a and s a and a spin- weight of +1 to a a and t a [22-24] . This allows us to 
express the results in terms of spin-weighted quantities which act as a check on the calculations and simplifies 
their appearance through the use of the edth operator which is also defined in the Appendix. In (4.4) and 
below, we exhibit only the terms of the solution we have calculated and omit the dots indicating further 
terms. 

The radial integrations are not unique, but lead to the introduction of a number of r— independent 
functions, C 1 s = C 1 a s a ,C 1 t = C x a t a , M, A, °V, and 1 V. The first four of these functions are related to the 
angular momentum, mass, and radiation. The remaining two are fixed by the reality conditions. Below we 
give the results of these radial integrations in the order given in Eqs. (3.2): 

=[^-> 1 £ + > 1 *-> 3 £K+ 

[- C ^ + ^\-^A] ba (4.5a) 
=^l^- 1 -C\ + ±(^-C\)] (4.56) 



S 3 a A 3 a = irsint 
A 3 ! - 0. 



(4.2a) 
(4.26) 



1 r 3 A 



•} 

3£ 



ir sin 9 .„ £ 
^{1 + — + — + •••.} 

t/2 r r 3 



(4.3) 



3A , 

-da + -^3»a 

-irsm9{s a + —t a + -^t a }. 



(4.4) 
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A\ = [^+i^ 1 S-^(M 1 S + ^C 1 s + 2(^ 1 S) 2 + 1 S^ 2l S]a a 

+ [~\+ l -M-^C\]b a (4.5c) 

V a = --\(- 1 V+-^°V)s a + {°V 3 A+-^r 1 V 1 T,)t a } (4.5d) 
r L r 2r J 2r^ J 

B 3 = 0(1) (4.5e) 

iV a = - 1 [| Vs a + ^f] (4.5/) 

iV 1 = - 1{°V- -[2 VA^ + ^ 2l V + ^( V^ 1 S + 2^ 1 V)]} (4.5g) 

Sl =^°V(f^-M) (4.5ft) 

^ 2 = + 1 V + ^(°V^ 1 S)]- 1 V^} (4.5i) 

The integration for iV a can introduce a function independent of r. However, such a term can be removed by 
a coordinate transformation [11,12]. Furthermore, N 1 will grow like r at null infinity unless we require that 

= = 0. 

The above solutions have been written with this requirement so that at most, °V can be a function of i alone. 

The fact that the constraints form a closed system shows that the propagation of these equations is 
consistent. This means that the propagation equations will determine the evolution of the arbitrary functions 
we have introduced, but there will be no further conditions. This argument is equivalent to the use of the 
Bianchi identities by Bondi and Sachs [11, 12] to obtain a similar result. 3 A, 1 S, and ^ are part of our initial 
data. In addition we have C 1 a , A4, A, °V, and 1 V. These are exactly the same quantities we would have had 
to introduce if we were to integrate the equations without the asymptotic expansion. 

From the evolution equation (3.3a), we obtain the following relations: 

°VA = "£ (4.6a) 



C s = -pB 1 + 2 B 2 + °V[-^(# 1 S-ls)+^M] (4.66) 

C t = -$ 2 B 1 + 1 V$A (4.6c) 

These latter equations are identified with the change in dipole aspect and, hence, are connected with the 
conservation of angular momentum. The equation for A 2 ^^ defines m which leads to an identity with Eq. 
(2.18). is the null part of the conformal tensor, ^4 in the Newman-Penrose notation [13]. Thus, only 

the equation for A 2 a V a is dynamical. Remembering that fiiV 1 = 0, the relevant relation can be written as 

-^(VM-V^) - (4.7) 

All that remains now are the equations to propagate A 3 a and £3°. This completes the solution of the 
complex equations. The further discussion of (4.7) will be delayed until after the reality conditions have 
been applied. 
The Reality Conditions. 
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It is only necessary to apply the reality conditions to the tetrad because through the solution of the 
field equations the connection is expressed in terms of the tetrad. We shall see that in applying the reality 
conditions, the arbitrary functions in the connection will be expressed in terms of those in the triad. The 
reality conditions on the tetrad are 

N = N, N* = N\ -iv = iv > 0, 

Ei* = -Si', v l = t 3 \ (4.8) 
From (3.1) we find that Si 1 already satisfies the reality condition. Writing v a = E 3 , we have 



E 3 = - ir sin 6{s a + ^Ef + ^t a }, 
v a = irsin6>{t a + J 1 Es a + ^L 3 Es a }. 

From this we find that 

v 1 = TfefcEiVEs* - rSin^fl-l^E]. 
The requirement —iv > gives us 



ir 



Now from Nv a — V a and the boundary condition that N should be 1 at null infinity, we find 



which implies N = 1 and 



Comparison with the Bondi-Sachs form of the metric [12,22] 

Ve 2b 

ds 2 = du 2 + 2e 2b dudr-r 2 h l:j (dx l -U l du)(dx3 -U j du), 



(4.9) 



2 sin(9Wl- ( 4 - 10 ) 



°V = -1, 1 V = - 1 E, ^ 3 A = 3 Z + ^ 1 E 1 T, 1 E. (4.11) 



(4.12) 



give us 

v = -ie 2b r 2 sinO, N = 1, iV 1 = — , 

r 

and for the principal spin coefficients 

p = - -A 3 a E 3 a = g-, 

v r r a 

v r A r 4 2 

A* = -^ 2 a« a = -^ + ^(M- 1 E 1 E), 
f 2r r" 1 

A = --A 2 a E 3 a = -- + ^(# 1 E-i 1 E). 
i/ r r z 2 

This shows that the Newman-Penrose asymptotic shear a = — E and 

-M + 1 E 1 E = ^ [V>° + V>° + # 2 ct° + ^ 2 ct° + 1 E 1 E + 1 E 1 E]. (4.13) 
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Note that the right hand side of (4.13) is real and comparison with (4.7) shows that it is just the negative 
of the Bondi-Sachs mass aspect. 

It is perhaps wothwhile to exhibit the real mass loss expicitly by writing the integral over a sphere at 
null infinity as (dS is the area element of the unit sphere) 



(4.14) 



M := — j [M + M - ^E - ^E] dS 
= /[X->(- 1 E 1 E + 1 E 1 E] 



(4.15) 



This definition of the mass aspect agrees with that of Bondi-Sachs and Newman-Penrose. Thus, (4.14) 
describes the mass loss from gravitational radiation. It is interesting to note that addition of the surface 
integral 

Iv a A 2 a i: 1 1 ded^ (4.16) 

is needed to assure the differentiability of the Hamiltonian. Together with the reality conditions, (4.16) 
defines the mass M. 

5. Logarithmic Behavior. 

In this section we shall consider the terms in 1/r which were omitted in the definition of (A 3 a ,Y, s a ) 
in the previous section. These terms lead to logarithmic behavior which comes in below the leading terms 
previously found. Chrusciel, MacCallum, and Singleton [15] have studied the logarithmic behavior within 
the Bondi-Sachs formalism [11, 12]. They introduce polyhomogeneous functions in the metric and then see 
what powers of (In r ) m /r n lead to consistency in the solution of the constraint and propagation equations. 
Here we follow a somewhat different approach. We start with a power series in 1/r in (A 3 2, E3 2 ) and examine 
the logarithmic terms which arise in the remaining terms. We then look at the propagation equations to see 
what logarithmic terms are introduced into (A 3 a , E3 ). Consistency means that these new terms should not 
interfere with what has been previously required. That is, a new calculation with these new terms should 
reproduce the already found behavior and add lower order logarithmic behavior. Apart from this difference 
in approach, our results are essentially in agreement with those of Chrusciel, MacCallum, and Singleton, but 
perhaps they are more perspicuous. Our approach is closer in spririt to that of Winicour [25]. However, in 
our results the coefficient of 1/r in A 3 a is found to be zero, whereas that does not appear to be the case 
in their analysis. For possible clarification, the details of the calculation leading to that result is given in 
Appendix 3. Below we sketch our final results and then discuss the conservation equations and the conformal 
tensor. 

From Gi,02, and Hi we find 



2 A 3 A 
A 3 a = -a a + -^b a + - T 



E 3 a 



-irsinO 



2S 



s a + —t a + -^t a + 



(3E - WA) 



(5.1) 



The later integration for V a requires that a possible term in 1/r in A 3 a be set equal to zero. Then the 
remaining integrations along the null generators give 



A 1 



cote 



2V2 

COt6» 



s 2 ^lnr 1 q X 
2\/2 r 2 r 2 * 

l x E + 2^ — 

r 



b a 



a a + 



(5.2a) 
(5.26) 
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A\ = 



~ + ^-^(f 2 A-^ 22 E) & a (5.2c) 



To the order considered, the solution for V a is the same as in the previous section. However, as noted above, 
it imposes the condition that A = 0. This result is contained in the above expressions. 

In the remaining variables, the logarithms appear in the order below the leading order given in the 
previous section. Except for N 1 , they have no important consequences. For N , we find a term in In r/r 2 . It 
then follows from the propagation equations that A 3 a develops a term in lnr/r 3 and T,s a a term in lnr/r 2 
which follows the term in X E . As a result, the conservation equation for mass is unchanged. This means 
that what appears as gravitational radiation at null infinity is unaffected by the inclusion of this logarithmic 
behavior. On the other hand, because A 1 a has a term in lnr/r 2 , the conservation equation for angular 
momentum will be changed. This represents another problem for angular momentum which is yet to be 
understood adequately. This is exhibited below by the failure of the conformal tensor to peel: 

r 



\E»4 


= R 2 pa eo p e 3 a 




= R 1 P ae 3 p e ° 




= i?Vei p e 2 ff 


*1 


= R 1 P aei p e 2 a 


*0 


= R 1 P aei p e 2 a 



(5.3) 



2 

„4 



lnr 2C a t - $ 2 A 



r* 1 

2 A lnr 3 A - a 

In the above, the \f n are the components of the conformal tensor, X S = — a° of Newman-Penrose, and a 
comes from the lowest order logarithm in A 3 a . 

The difference between our results and those of Chrusciel, MacCallum, and Singleton [15] comes from 
the different question which is asked. We ask for the logarithmic behavior which is forced on us by adding 
the additional terms in 1/r in both A 3 a and S3 while they ask for the most general logarithmic behavior 
they can introduce into the metric which is self-consistent. As a result, they find that the metric can have 
terms in whose coefficients are independent of u. However, the main physical conclusions arc the 

same. The definition of Bondi mass and the radiation of gravitational energy remains unchanged from the 
results without logarithmic behavior. Furthermore, considerations about angular momentum are affected by 
these new terms. On the other hand, the current results are in complete agreement with those of Novak and 
Goldberg [16, 17] who showed that the existence of null infinity was consistent with the logarithmic behavior 
found here. 



13 



6. Conclusions. 

One of the main points of this paper has been to see whether the new variables introduced by Abhay 
Ashtekar [1-4] are useful in classical problems. Indeed, use of the self-dual connection and the densitized 
triad (A- A a ,T,_^ a ) in the canonical formallism leads to a set of equations which are intermediate between 
the Bondi-Sachs [11, 12] and Newman-Penrose [13] equations in the vicinity of future null infinity, I + . In 
Bondi-Sachs, the calculation begins with specification of the metric on a two-surface foliation of an outgoing 
asymptotic null cone; in Newman-Penrose, with the specification of a component of the conformal tensor, 
^oN 2 Cijk m ^i l ^3 : '^i k ^3 rn ] and in the present calculation we give as our initial data (A 3 2,T,s 2 ) which is 

only half of the two-surface metric and that part of the connection which is related to ^o- We have not 
evaluated whether our calculation is the most efficient. Both Newman-Penrose and we work with first order 
equations. Eventually they make use of the "metric" equations to determine the tetrad whereas they are 
part of our canonical equations, but we must then compute the conformal tensor. 

The one advantage of the present approach is that it is derived from a Lagrangian and a canonical 
formalism which may yet be useful for quantum gravity. It also makes the study of the logarithmic behavior 
in the vicinity of null infinity somewhat more imperative and somewhat easier. Apart from the point 
mentioned earlier and elaborated on in Appendix 3, the important conclusions we have arrived at are not 
significantly different from those of Chrusciel et al [15], if less complete. We hope that this difference can 
be resolved in the near future. As noted in the previous section, we ask different questions. The important 
results here and there are that one can have an asymptotically Minkowskian metric, with a future null 
infinity, and a mass and radiation of gravitational energy which is well defined. That is, the logarithmic 
terms fall off faster than those terms which define the mass and radiation of gravitational energy. The same 
is not true for angular momentum. But, that concept is not sufficiently clear even in the absence of the 
logarithmic terms, although, in that case, there is an expression for angular momentum which transforms 
properly under the super translations as well as the Lorentz transformations [25]. 

There is, in the present approach, the additional need to apply the reality conditions. However, note 
that there are no propagation equations for v l . Therefore, once v l is identified with S 3 J , it follows that its 
propagation is also specified. Given that £i J is real and independent of time, the metric will propagate as 
real. That guarantees that all the reality conditions will be fulfilled. 

Our original hope was that on the null surface we could carry out a reduced phase space quantization of 
general relativity. While the identification of the dynamical degrees of freedom for the gravitational field is 
easy, either one has to express all the remaining variables in terms of these degrees of freedom or construct 
the Dirac brackets. At this time, it appears to be very difficult to carry out that task. 
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Appendix 1. Minkowski Space Tetrad and Connection. 



In Minkowski space the metric takes the form 

ds 2 = dt(dt + 2dr) - r 2 (d6 2 + sin 2 6 d<j> 2 ) 
so that the four one-forms and tetrad vectors are 

e° = dt, e = 8t - \dr, 

e 1 = (dr + \dt), ei = d r 

02 = 7f r (^ - isinOdcP), e 2 = j=(d e + ^0^) ( AL1 ) 

03 = 72 r ( dd + ^in^#), e 3 = ^ 2 (d e - ^dj,) 

Therefore, the triad densities and the self-dual connection one-forms are 

Ei* = —ir 2 sin# 5 i i; A 1 i ^cos96 l 3 ; 

v a =ir sine t a ; A 2 a = -\b a ; (A1.2) 

S 3 a =-irsin6s a ; A 3 a = -a a ; 

a a = -4=(1, isind); 
V2 

b a = ~isin9). 



(AIM) 



d 


i d ^ 


86 + 


sin# <9<^/ 


(- 


i <9 


\de 


sin 9 deb 



The duality relations between s a , t a and a a ,b a , respectively, are unchanged by a phase change e~ lc ^s a , e l ^t a and e l ^a a , e^^t 
We assign a spin-weight of — s to those quatities which transform as e~ ls ^ and a spin-weight of +s to those 
which transform as e ls ^ under this possible change of phase. Keeping track of the spin-weight is a help in 
controlling the calculations [22-24]. 

This also leads to the introduction of spin-weight raising and lowering operators $ and <f . If r\ has 
spin- weight s, then 

h--= sin* 0(1; + ^^(sin-^r?) {AIM) 

and 

h--= -^0(^ + ^^Or,). (AIM) 

The action of these operators on the spin-weighted spherical harmonics can be found in Appendix A of [22] 
and further details can be found in [23,26]. 

Appendix 2. Properties of Z a b 

In Eq. (3.2e), we defined (a,b = 2,3) 

Z\ := 5\ - A, (,42.1a) 
ft - §£j 

It is easy to see that the determinent of /% is zero and there are two eigenvalues (1,0) with the respective 
eigenvectors £3° and e ab Af and eigenforms A 3 a and e a fcS3 b where e a b is the two-dimensional Levi-Civita 
tensor. The same is true of Z a b with the eigenvalues (0,1). Therefore, Z a b is a projection operator. In the 
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1 jr expansion of the variables, the zeroth order term will select eigenforms and eigenvectors of spinwcight 
— 1 and + 1, respectively. More specifically, we have 



1 / 1 —is'm9 

2 \i/sm9 1 



(A2.2) 



so that 

a Z\s b = a a a Z\ = 

042.3 

°Z\t h = t a , b a Q Z\ = b b 

The vectors and covectors s a 1 t a ,a ai and b a are the spin- weighted basis vectors and covectors introduced in 
Appendix 1. 

Appendix 3. The Solution for V a . 

We make no restrictions on E 3 2 or A 3 2 except that they should start with their Minkowski space values 
at null infinity and have a 1/r expansion in the vicinity of null infinity. In addition there are the coordinate 
and tetrad conditions which are described at the beginning of Section 3. We could proceed as in Section 4 
specifying S3 2 and A 3 2, but equivalently we can write 

S 3 a = -irsin9{s a + + 4 2 S a } 

! r ! r (^.1) 

A 3 a = -{a a + - x A a + -j 2 A a } 
r r z 

The restrictions which are a result of the equations Gi,02, and Hi lead to the same result. G2 tells us that 
A 3 1 — 0, then Gi and Hi give us 

Z 3 a A 3 a = ir sm9 

(A3.2) 

A 3 aA X 3 a = 0. 
Substituting (A3.1) into these equations yields the result 

E 3 a = -irsin6{s a + -^f + -^( 2 Zt a - I^Ss"} 

' r2 (A3.3) 

A 3 a = -{a a + - l Ab a + -( 2 Ab a - ^A^a a } 

From this point on, one just puts this result into the succeeding equations and proceeds as before looking for 
the solutions. Now some logarithmic terms appear, but below the first couple of terms. However, the problem 
arises in the equations for V a ,x a — 0, which do not depend on these solutions. We have (V a := Nv a ): 

(y^x 1 ),! - ^Z a b (V b ^ 1 1 ) = 0. (.43.4) 



Z\ := S a t 



^3 a A 3 b 
Z 3 C A 3 C 



Write 



:= t a b b - 1 [s a b b x A + t a b b 1 Y 1 ] - ^[{s a a b +t a b b ) 1 A 1 ^ + 2 Y]. 

V 11 ^ 1 = -irsm8[°V a + — + ^J] (A3. 5) 



r 

and substitute into (A3. 4). We obtain 

V - ^ = [rt t -J(sV^ + t a b b 1 E)-^[(s a a b + t a b b ) 1 A 1 E + 2 T,] 

lyo 2ya 



(A3.6) 



x [V + — + -V] 
16 



which gives the relations 

°V a - °V b b b t a , (V a - < V 6 a& 1 £)t a - VfyU s a - 0. 

Thus, 

°V a a a = Vbt = (°V a 6 a )U = 0. (3.7) 

The other expansion terms are given in terms of known quantities except for °V a b a and 1 V°a a which remains 
arbitrary. In order that v a = £ 3 a when the reality conditions are imposed, we necessarily choose X A = as 
noted in Section 5. Then we have 

2/ 34 (^3.8) 

-4 3 a = - a a + —b a + -^rba- 
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